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THE SECULAR DISPLACEMENT OF THE ORBIT 
OF A SATELLITE. 



BY PROFESSOR ASAPH HALL. 

(1). In his Mecanique Celeste, Tome II, p. 373, Laplace has considered 
the cause that holds the Rings of Saturn nearly in the plane of the equator 
of the planet, and he has explained the reason that led him. to predict 
the rotation of Saturn and its Rings before Herschel had determined these 
rotations by observation. In this investigation Laplace employs the equa- 
tions that determine the rotation of a solid body, the disturbing forces being 
those that arise from the action of the excess of matter around the equator 
of the planet and the attraction of a distant body. In the same chapter it 
is also shown why a satellite will be held nearly in the equator of its pri- 
mary when the figure of the planet is elliptical. These results have been 
used by Professor J. C. Adams to explain why the orbits of the satellites of 
Mars are at the present time nearly coincident with its equator. M. F. 
Tisserand has discussed the same question by means of formulse derived by 
him in a discussion of the secular perturbations of the orbit of Japetus, the 
outer satellite of Saturn. The investigations of these astronomers render it 
probable that the figure of Mars is slightly elliptical, and that the orbits of 
its satellites will always remain nearly coincident with the equat(U' of the 
planet. 

The question of the perturbations of satellites has been treated by La- 
place in the fourth volume of his Me'canique Ce'leste. There Laplace 
determines the perturbations of the radius vestor, and of the longitude and 
latitude of the satellite ; and for astronomical purposes these coordinates 
seem to be as good as any that can be chosen. If we wish to know, how- 
ever, the effect of the perturbation on any element of the orbit, it is more 
convenient to use the formulae of Lagrange. In the unfinished essay on the 
Saturnian System left by Bessel he has given these formulse, and also the 
various forms of the perturbative, or potential function that come into use 
in this complicated and interesting system. M. C. Souillart has lately pub- 
lished an elaborate memoir on the satellites of Jupiter in which he verifies 
and extends the investigations of Laplace referred to above. 

The forces that disturb the motion of a satellite are, (1), the action of the 
protuberant matter around the equator of the primary planet; (2), the at- 
traction of another satellite of this planet; (3), the attraction of a body out- 
side this system, like the sun or another planet; and in the case of Saturn 
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the action of its Ring also produces a disturbance in the motion of its satel- 
lites. Since the perturbations are small they may be computed separately 
and then added together ; I consider here only those which are produced by 
the protuberant matter around the equator of the planet, and by the action 
of the sun. 

Denote by d the longitude of the ascending node of the orbit of the satel- 
lite on the orbit of its primary planet, and by i the inclination of the satel- 
lite's orbit to the same plane. Let g be the angle between the line of apsides 
of the orbit of the satellite and its line of nodes ; then if a be the semi-ma- 
jor axis of this orbit, e the eccentricity, and if n be the mean sidereal motion 
of the satellite, we shall have by the formulse of Lagrange; 

<Ji a« /dR\ an cos i /dlt\ 

'(U~siQii/{l—e^)' \~W) sinVll—e") ' ^W) 

dl ^ m_ /dE\^ ^' 

dt sinij/(l — e^) ' \ di /' 

i2 denotes the perturbative function used by Laplace. 

We have now to expand this function for the two disturbing forces in 
elements of the orbit of the satellite, substitute the partial diiferential coef- 
ficients of i? in equations (1), and then integrate them. 

Assuming the satellite to be a material particle whose mass may be neg- 
lected, let X, y, z, be its rectangular coordinates referred to the centre of 
gravity of the planet, which is supposed to be fixed, and let 

Let X, Y, Z, be the coordinates of the centre of the sun, S its mass, and 

also let 

D = /(Z^ + Y^ + Z% 

If m be the mass of the primary planet and — ■\- Y the sum of the 

particles of the planet divided respectively by their distances from the sat- 
ellite, then the value of the perturbative function will be, 

7? =. S^Xx^Yy+Zz) S y 

D^ [j^X-xY^{Y-yY+{Z-zfY 

The first two terms give the perturbations produced by the sun, and the 
last term those caused by the non spherical figure of the primary. We have 
now to develop this value of R and find those terms that produce the secu- 
lar part of the perturbations. The denominator of the second term is 
D [lyT'^ _ ^{Xx+Yy +Zz)\l 

and expanding the reciprocal of this we get for the second term of JS, 
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__ ^ -L 'S-^' _ S-(Xx+ Yy+Zz) _ SS.r* 

S8.{Xx+ Yy+Zzf , SS.r\Xx+ Yy+Zz) „ 

The first term of this expausion may be omitted since it does not contain 
the elements of the orbit of the satellite and will disappear in the differen- 
tiation, and the third term cancels the first term in H. Again D is very 
great with respect to r, and hence the approximate value of the perturbative 
function that depends on the action of the sun becomes 

P _S.r' 38.iXx+ Yy+Zzf .„. 

^ ~ 25"' 25^ ^' 

If a great circle be drawn from the position of the satellite in its orbit to 
the sun, and we call the arc joining these bodies /, we shall have, 

. Xx+ Yy+Zz 

and therefore 

^ ^^.^_ {Xx+Yy +Zzr 

If we denote by u and U the angular distances of the satellite and the 

sun from the node, the spherical triangle between the node, the satellite and 

the sun gives 

cos / ^ cos M cos Z7 -|- sin M sin ZJcos i, 

and hence 

cos/^ ^ cosM^cos U^+ sinw^sin f/^cost^ 

+ 2 cos M sin u cos Z7sin ?7cos i. 
Changing the powers into cosines of the multiple arcs we have 

cos/^ = l^.p — sini^ + 2 cos 2m cos 2 C7 + 2 sin 2m sin 2 ZJcos i 
+ (cos 2m+cos 2 Z7— cos 2m cos 2 ?7)sin *'^]. 
As we wish to get only the secular part of this expression we have by 
omitting the periodical terms depending on thi; position of the sun, 

^,.cos/^ = -rro--^ (2 — sini^).r''+r^cos2M.sin*'' \. 
Ir 4:lr [ ) 

If now he a, function of the radius vector and true anomaly, r and v, 
in an ellipse, and if we denote by if and U the mean and excentric anoma- 
lies, then the non-periodical part of is given by the definite integral, 

1 r0.dM = ,^— 1 — ^ . r> .^.dv=^ f>. ^ . dE, 

Since from the equation of the ellipse we have 

r := a(l — e cos E), 
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the last integral gives for the non periodical part of r", 
^ C\l ~e cos E)^.dE= a\{l+-^^e'). 

Again since u = v-\-g we have to find the non periodical parts of 

= r^cos2»; and <? = r^sin2«. 

T he equation of the ellipse gives 

sin ^. -1/(1— e^) cos^— e 

sin V = ■ ^ ,-, : cos « = :- =, 

1— e cos ^ 1 — e cos E 

and hence 

r^ cos 2v = a\ [(cos ^— e)^ — (1 — e") . sin E^^, 

r^sm2v ^ 2a\smE .{cosE— e).i/{l — ^). 
Substituting these values of # and integrating we have 

= fa^e^- and = 0. 

If flfl and Bq be the semi-major axis and the excentricity of the orbit of 
the planet we shall have for the non periodical part of -=yg, 

#==j_ r II = I 

27ral -' o(l— eoCosjB)2 a^^{l—el)i 
From these results we find for the secular part of B^, the value 

The term cos-=^ in cos /^ produces no constant terms of an order lower 

than e^, and therefore it may be neglected. 

(2). In the case of an ellipsoid of revolution Laplace gives a very sim- 
ple form to the potential V of the disturbing force. If we denote by p the 
ellipticity of the spheroid, by cp the ratio of centrifugal force to gravity at 
the equator of the planet, by B the radius of this equator, and by /j. the sine 
of the declination of the satellite with respect to the same equator, we shall 
have; Mec. Cel., Tome II, p. 103, 



V = 



Taking the mass of the planet and its equatorial radius for the units of 
mass and distance, and putting Q' = h<p — p, we have 

Let </) be the distance from the node of the orbit of the satellite to its in- 
tersection with the equator of the planet, and let y be the inclination of its 
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orbit to this equator, then from the right angled spherical triangle we have 

/x := smy sin {u — <p) 
and 

/z^ = J sin y^ . [1 — cos 2{u — ipj]. 

The part of V which gives the secular terms is therefore 

Let A be the inclination of the equator of the planet to the orbit of the 
planet, and let d^ be the longitude of the node of the equator on this orbit; 
the spherical triangle formed by the three nodes gives, 

cos ;' = cos (6 1 — d) sin J. sin i + cos A cos i 
sin ysin (p = sin {6^ — 6) sin A 
sin ;' cos (p = cos {6^ — 6) sin A cos i — cos A sin i. 

These equations give the means of expressing sin ■f in terms of d and i, 
and hence i?2 becomes a function of the elements of the satellite's orbit. 
The values of J. and d-^, which determine the position of the equator of the 
planet on the plane of its orbit, are in fact subject to small variations simi- 
lar to precession and nutation, but here these quantities are supposed to be 
constant. 

The excentricities of the orbits of the satellites in our solar system are 
generally very small, and if we neglect the squares of these excentricities, 
and omit the constant part in i?]^, equations (1), (3) and (4) give, if we put 



q = 



8ag(l-eg)f 
an fdR \ 

\di) 



dt sin i \dd / 

dd an 

dt sin i 



or we may write 
Hence we find 



V di 
The value of cos y gives 



0' 
B =Q. sint^— ^ . sin y^; 

B = ^e. sin i^ -(- Jc' . sin y^. 

/dB\ . . . , fd.cosr\ 

1—^1 =c.smicosi — c'.cos;'.! — ' 1. 
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( ' 'J = sm(^j — ^) sin J. sin i 

i ' ' ] = cos(^] — ^)sin^cosz — cos J. sin t. 
The other equations of the triangle enable us to change these values to 
I — -j^ 1 = sin ;- sn: </> sin i 

I — ;p-^) = sin;- cos ^. 

We have therefore finally 

di ,, ... 

-j- = — k. cos J- sin Y sin ^ 



dt 
dd 



(5) 

, = — k cos i -4- h'. - 
dt ^ 



cosj- 


sm j-i 


BOS^ 




sin* 


y 


¥■■ 


:=(£: 


«.2 



where h and ¥ have the values, 

4a3(l— e§)l 

These are the formulae given by Laplace, Mec. Cel., Tome IV, page 182, 
The quantity h can be computed accurately for all of the principal planets ; 
but ¥ is not so well known since it depends on the ellipticity of the figure 
of the planet and on its law of density. If the planet has other satellites 
the secular perturbations of these should be comprised in the value of ¥, 
and in the case of Saturn ¥ would include also the secular action of the 
Ring- 

The Italian astronomer Plana criticised equations (5), and considered 
th(im erroneous; but in the Conn, des Terns for 1829, p. 245, Laplace has 
verified his first proof of them, and has given also the simple derivation 
from the formulae of Lagrange which has been followed above. 

(3). In his investigation of the motion of the orbit of Japetus Laplace 
has introduced an auxiliary fixed plane passing through the line of nodes 
and placed between the equator of the planet and its orbit, the angle which 
this plane makes with the equator of the planet depending on the ratio of 
k to ¥. 

M. Tisserand has employed a property of the perturbative function which 

leads very simply to a knowledge of the curve described on the heavens by 

the pole of the orbit of the satellite. If £ be the mean longitude at the epoch, 

and ;r be the longitude of the inferior apsis, the complete differential of B is 

JO dR , , dB J , dB -j^ , dB -j^idB j-,dB ,n 
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but for the secular perturbations we have 



and hence 



da = 0, and —r- = 0, 

as 



dR=^.de+^ .dn+^-^ .d.i+ ^ .dd. 
de uTc di da 



Now 


we have 


by 


the theory 


of Lagrange, 










dR 

de 


(e, tt) . 


S+(.^). 


dd 

dt 








dR 

dlT 


(JT, e) 


de 
■ dt 










dR 
di 


a, d) 


dd 
' dt 










dR 

dd " 


{e,i) 


.|me) 


de 
' dt' 



If we substitute these values of the partial derivatives of R in the second 
value of dR and notice the well known relations 

(e, 7r) = — {n, e), &e., 
we shall have 

dR==0, 

and hence 

R = constant. 

From the value of R that we have found we have, 

k sin i^ -|- k' sin y^ = constant. (6) 

If a satellite therefore be acted on by the two disturbing forces that we 
have considered the pole of its orbit will describe on the heavens a spherical 
ellipse. This equation shows also that if the disturbing force of the sun be 
zero, the orbit of the satellite will have a constant inclination to the equator 
of the planet, and if k' vanish, it will have a constant inclination to the or- 
bit of the planet. The real position of the orbit of the satellite will depend 
on the ratio of the quantities k and k', and for this ratio we have 

or, since the mass of the planet is the unit of mass, and 
1 = a^n"^, (S'= a^nl ; 

where the unit of distance is the equatorial radius of the planet. The un- 
certainty of this ratio lies in our ignorance of the factor (p — J^), For the 
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planets Mars, Jupiter and Saturn the quantity f is known with a tolerable 
degree of accuracy, since it depends on the time of rotation of the planet 
about its axis ; but the value of p must be inferred from the probable con- 
formation of our Earth, or from some other analogy. In the case of the 
exterior planets of our system the ratio k'-^lt is probably large for all the 
satellites, with one exception, and for this reason these satellites are found 
nearly in the equators of the planets. This ratio, it will be seen, varies in- 
versely as the fifth power of the mean distance of the satellite, and it is on 
account of this fact that Japetus, the outer satellite of the Saturnian sys- 
tem, presents an exception to the general rule, its orbit being inclined 
nearly fourteen degrees to the equator of Saturn, while the next interior 
satellites, Hyperion and Titan do not depart more than one or two degrees 
from the equator. In the case of Japetus Laplace found, 

V 
-f = 0.4219, 

which is too large on account of the uncertain data used in the calculation. 
M. Tisserand finds from Bessel's value of the mass of the Ring and ellipti- 
city of the planet, 

^ = 0.2570; 

and this ratio is still uncertain since the mass of the Ring is not well known, 
and the masses of all the satellites of Saturn are unknown. The better way 
probably would be to determine k' by observations made at two distant 
epochs which would enable us to fix the position of the node at those epochs 
and the values of the inclination, and by comparison of these values we 
should have k'. In order to distinguish in the case of Saturn the different 
parts of k' it will be necessary to determine accurately the motions of the 
interior satellites, since from these motions the mass of the Ring and the 
effect of the ellipticity of the planet can be found. 

For the application of equations (5) to Japetus we need the values of the 
inclination and longitude of the node of the orbits of Saturn and Japetus 
and of the equator of Saturn on the ecliptic. These values are for 1880.0 

N. J. 



Saturn 


112° 


36'.5 


2° 


29'.6 


Japetus 


142 


45.0 


18 


31.5 


Equator of Saturn 


167 


55.2 


28 


10.3 



The solutions of the three spherical triangles formed by the nodes of 
these great circles give 
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i = 16° 


25'.0 


r = i3 


41.0 


^ = 53 


39.0 


A = 26 


49.6. 



The periodic times of Saturn and Japetus are 

To= 10759.2198 days 
T = 79.32936 " 
The value of h may be written, 



(1-62)1 

For Saturn e^ = 0.0559956; and if we take the Julian year of 365.25 
days for the unit of time we have for Japetus n = 5967071".2. The val- 
ue of k is 

k = 0.0000409651 . n. 

Equations (5) give for the annual variations of i and d, 

S = ~ 45".250.-^ 
k 

AO = — 234".476+117".827 .—. 

k 

With Tisserand's value oik' ^ k the annual motion of the node on the 
orbit of Saturn will be — 204'M94. 

For the outer satellite of Mars the values of the nodes and inclinations on 
the ecliptic for 1880.0 are: 

N. J. 

Mars 
Deimos 
Equator of Mars 

With regard to these quantities it must be noticed that the position of the 
equator of Mars is uncertain; but from the data given we find, 

i = 24° 20'.0 

r = 4 5 .6 

(J) = 312 38 .5 

JL = 27 15 .9 

The periodic times of Mars and Deimos are 

ro= 686.97965 days 
r = 1.2624350 « 
and in the orbit of Mars e^ = 0.09328975 ; and for the satellite 
log n = 8.5739862. 



48° 


37'.9 


1° 


51'.0 


85 


37.4 


25 


47.2 


79 


89.0 


28 


60.0 
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The value of k Is 

h = 0.0000025661722. n. 

and equations (5) give for the annual motions of i and d on the orbit of 

Mars: 

Ji = + 60".396 . ^ 

Ad= — 876".735 + 112".630 . ^. 

The node would have therefore a retrograde motion by which it would 
complete a revolution on the heavens in 1478 years were it not for the in- 
fluence of the second term. In the case of Mars the value of the ratio It -irk 
is unknown, but probably it is large enough to destroy the first term in the 
value of AO. 

(4). To find the equation of a spherical ellipse, let 2a be the part of a 
great circle forming the major axis of the ellipse, and let 2e be the distance 
between the foci. Denote by p and p' the arcs drawn from a point on the 
ellipse to the foci, and we shall have the condition, 

p -\- p' = 2a. 

If we take the centre of the ellipse as the origin of the polar spherical 
coordinates S and F, S being the mdius vector and V its angle with the 
major axis, we have, 

cos p = cos ;8'cos £ — siu (8 sin e cos V 
cosp' = cos iS cose + sin/Ssin ecos V. 

Adding and subtracting these equations, and reducing, we find 

cos Scos e 



cos Up'—p) = 



cos a 



. 1 / / * — sin /S sin e cos V 

sin i(p' — p) = ; ; 

^^'^ ^^ sin a ' 

and squaring and adding these we have for the equation of the ellipse, 

1 = 22?_^ . cos ,S'+?i5-^, . sin 8' COS V. (8) 

cos a* sma^ ' 

The form of our perturbative function is the same as that of equation (8) 
since it may be written, 

B = keoa'P + k' cosy' — c. 

The spherical ellipse will be symmetrical with respect to the great circle 
joining the pole of the orbit of the planet to the pole of its equator, and 
therefore the centre of the ellipse will be on this circle and between these 
poles. To find this centre M. Tisserand proceeds as follows ; let vj and rj' 
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be the distances of the centre from the pole of the orbit of the planet and 
from the pole of the equator. Then 

7] -{■ yj' = A. 

We shall have from the spherical triangles 

cos i = cos 7j' cos 8 — sin yj' sin 8 cos V 

cos ;' := cos ^ cos /S -[- sin ;y sin 8 cos V. 

The equation of the ellipse will be 

A.(cos ■^' cos 8 — sin tj' sin -S cos Vf -f i!;'.(cos;y cos /S'+sin r) sin Scos Vf = c. 

By equation (8) the coefficient of sin /S cos iS' cos Fmust be zero. Hence 

— 2fc.cos 3y' sin tj' + 2^'.cos r^sinvj = 0, 
or 

Lain 2{A — jj) = h'.sm tq, 
and 

n J.sin 2 J. 

*"°S^'^ = F+I^3^23- 

The angle t] is the same as the angle d used by Laplace, Me'c, Ce'leste 
Tome IV, p. 178. For two satellites of the same planet the ratio 1<! -^ h 
will vary inversely as the fifth power of their mean distances. If we assume 
for Japetus this ratio to be 

4^ = 0.2570, 

A! 

its value for Ehea, the fifth satellite of Saturn, will be 

Hence sy = 22".9 ; or, for an observer on our Earth Ehea will moye 
sensibly in the equator of Saturn. For Titan, the large satellite next be- 
yond Rhea, we have 

T - ^'^^^^ V22J45J • 

Hence ;y = 1520".5. This is smaller than the value given by observa- 
tion, and indicates that the assumed ratio 

I' 



^ = 0.2570 



is too great. 



